Nonlinear Propagation of Incoherent Photons in a Radiation Background 
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The nonlinear propagation of intense incoherent photons in a photon gas is considered. The 
photon-photon interactions are governed by a pair of equations comprising a wave-kinetic equation 
for the incoherent photons in the presence of the slowly varying energy density perturbations of 
sound-like waves, and an equation for the latter waves in a background where the photon coupling 
is caused by quantum electrodynamical effects. The coupled equations are used to derive a dis- 
persion relation, which admits new classes of modulational instabilities of incoherent photons. The 
present instabilities can lead to fragmentation of broadband short photon pulses in astrophysical 
and laboratory settings. 

PACS numbers: 12.20.Ds, 95.30.Cq 



The advent of quantum electrodynamics (QED) is one 
of the major scientific achievements of the 20 th century. 
It is an experimentally well confirmed theory, and it has 
predicted a number of phenomena which were not ex- 
pected in previous studies, e.g., the Casimir effect. More- 
over, as opposed to the classical theory of Maxwell, elec- 
tromagnetic waves can according to QED interact in the 
absence of a material mediator. Due to the possibility of 
exchanging virtual electron-positron pairs, there is the 
intriguing effect of photon-photon scattering, first dis- 
cussed even before the interaction between light and mat- 
ter was well understood 0, , and later derived within 
QED by Schwinger Q . The fact that strong electromag- 
netic fields can interact in vacuum, opens up for very 
interesting applications in extreme astrophysical environ- 
ments, such as magnetars |4L where photon splitting and 
lensing may take place [El* HI • The prospect of di- 
rect detection of the effect has also been discussed in the 
literature, and a number of suggestions for experimen- 
tal setups have been given, involving second harmonic 
generation self-focusing |10|. nonlinear wave mixing 
in cavities 0, 0] and waveguide propagation 0] , re- 
spectively. Apart from the astrophysical regions where 
high fields exist, the rapid development of ultra- high laser 
fields [l4l llil and related laser-plasma techniques (lt| . 
gives hope that the critical field strengths will be pro- 
duced in laboratories. 

We shall below investigate the effect of incoherence of 
high-frequency photons propagating on a radiation fluid 
background. This will extend the work presented in Ref. 
[T?| . and investigated numerically in Ref. [l^, to the 
random phase regime. For this purpose, a wave kinetic 
theory for high-frequency photons coupled to an acous- 
tic wave equation for a radiation fluid is presented. It is 
shown that modulational instabilities are inherent in the 
system of equations. Moreover, in the limit of the slow 
time-variation approximation, we obtain a Vlasov equa- 
tion with self-interaction for the high-frequency photons. 

We consider an incoherent non-thermal high-frequency 



spectrum of photons. As will be shown, this spectrum 
will be able to interact with low-frequency acoustic-like 
perturbations. The high-frequency part is treated by 
means of a wave kinetic description, whereas the low- 
frequency part is described by an acoustic wave equation 
with a driver |17| which follows from a radiation fluid 
description. Let Nk(t,r) denote the high frequency pho- 
ton distribution function, normalised such that the cor- 
responding number density is given by n = J N k d 3 k. 
Then N k will satisfy the wave kinetic equation [ljj 



dN k 
dt 
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where v g = du k /dk is the group velocity, and 
co k = ck (1 - , 
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where 6> is the radiation fluid density |17tl2C|| . and A = 8k 
or 14k, depending on the polarisation state of the photon. 
Here k = 2a 2 ti 3 /45mjc 5 m 1.63 x 1CT 30 ms 2 /kg, a is the 
fine-structure constant, 27r?i the Planck constant, m e the 
electron mass, and c the velocity of light in vacuum. The 
dispersion relation (J2J) is valid as long as there is no pair 
creation and the field strength is smaller than the QED 
critical field, i.e., 



lu < m e c 2 /h and \E\ < E t 



crit 



m e c 2 /e\ c 



(3) 



respectively. Here e is the elementary charge, A c is the 
Compton wave length, and E CT it — 10 18 V/m. 

The high-frequency photons drive low-frequency 
acoustic perturbations according to |17| 



dt 2 



2X(OQ 



dt 2 



c 2 V 2 



hui k N k d 3 k 
(4) 

where the constant §q is the background radiation fluid 
energy density. This hybrid description, where the 
high-frequency part is treated kinctically, and the low- 
frequency part is described within a fluid theory, ap- 
plies when the mean-free path between photon-photon 



FIG. 1: The transverse instability for the mono-energetic case. 
r/Kc plotted as a function of vt / c, as given in Eq. JHJ . 



FIG. 2: (F '/ 'Kvt) 2 , according to Eq. plotted as a function 
of qo = cos 80 in the mono-energetic case. 



collisions is shorter than the wavelengths of the low- 
frequency perturbations. We note that the intensity 
I k = hu> k N k /eo satisfies Eq. and is normalised such 
that (l-El 2 ) = J Ikd 3 k, where E is the high- frequency 
electric field strength, and eo is the dielectric constant 
of vacuum. The equations presented here resemble the 
photon-electron system in Ref. |2l|, where the interac- 
tion between random phase photons and sound waves in 
an electron-positron plasma has been investigated. 

Next we consider a small low-frequency long wave- 
length perturbation of a homogeneous background spec- 
trum, i.e. N k = N kQ + N kl exp[i(Kz - Qt)], N kl < N kQ , 
and S = S\ e~xp[i(Kz — fit)] and linearise our equations. 
We thus obtain (using expression @ for tv k , and intro- 
ducing k = k/k) 



N kl 



2\kS l 



Kc 



ft - Kck ■ z 



dNko 
dk 



(5a) 



and 



Si 



2XcM Q f2 2 + K 2 c 2 
3 !l 2 -K 2 c 2 /3 



kN kl d 3 k, (5b) 



which, when combined, give the nonlinear dispersion re- 
lation 



2„2 



1 = - 



(j,K n 2 +K 2 c 

~n 2 ~K 2 c 2 /3j n-Kck-z" dk 



, (fi) 

where we have introduced the constant \i = ^\ 2 c 2 Ti<oq. 



As was found in Ref. 21, 22} , we may have growth for a 
large class of background distributions N k0 . 

(a) For a mono-energetic high frequency background, 
we have N k o — noS(k — fco). The nonlinear dispersion 



relation © then reduces to 

(n 2 - K 2 c 2 /3)(Q - Kccos8 ) 2 = ^n k K x 
x{H 2 + K 2 c 2 )[Kc + (2/2- 3ifccos^ )cos 6> ],(7) 

where we have introduced cos6>o = ko ■ z. This mono- 
energetic background has a transverse instability when 
8q = 7r/2, with the growth rate 



r =7E Kc 



-1 1/2 



/Vr_ 
V C 



14 



i-'T 1 I 



(8) 

where F = — ifl, and v T = /ingfcgc, and vt is a character- 
istic speed of the system. The expression in the square 
bracket is positive definite. In Fig. ^ we have plotted 
F I Kc as a function of vt / c. 

In fact, under most circumstances, vt <C c. Using the 
expression ((7J), we then have two branches. The branch 
corresponding to Q ~ Kc/^/3 is always stable for small 
vt, while for the branch corresponding to Q ~ Kc cos 9q 
we obtain the growth rate 



r = Kv T 



1 — cos 6*o 
1 — 3 cos 9q ' 



(9) 



which is consistent with (JSJ) in the limit 8q — > tt/2. In 
Fig. [21 the behaviour of the growth rate © is depicted. 

(b) The high-frequency photons have generally a 
spread in momentum space. For simplicity we here 
choose the background intensity distribution as a shifted 
Gaussian, i.e. 



IkO = 



-fa 



7T 3 / 2 k 3 w 



exp 



(fe - fco) 2 



k 2 



(10) 
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where — (l-^ol 2 ) is the (constant) background inten- 
sity and kw is the width of the distribution around k . 
Then the dispersion relation is 



1 = 



1 



k 5 w rj 2 - 1/3 



x exp 



k(kg cos#o — kcos( 



rj — cos ( 



k 2 



d 3 k, 



(11) 



where b 2 = (4/97r 3 / 2 )A 2 e <f ^o exp(-ftg/*£r) and V = 
Q/Kc. 

Assuming that the deviation of fco from the z-axis is 
small, and that 5 = ko/kyy <C 1, we can integrate Eq. 
Q11JI. keeping terms linear in 8, to obtain 



-irb 
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rj — 45r) cos 9q (2 arctanh?7 — in) 
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for < r\ < 1. Thus, we see that the non-zero width of 
the distribution complicates the characteristic behaviour 
of the dispersion relation by a considerable amount. It is 
clear though, that the width will introduce a reduction 
of the growth rate, as compared to the mono-energetic 
case. 

We may also look at the case when the time- 
dependence is weak, i.e. d 2 S '/dt 2 -C c 2 V 2 <f, such that 
Eq. J3J yields 



(13) 



S = 2X^0 J huj k N k d 6 k. 

Upon using Eq. (0 in ©, we find 

-uk^- [ k'N k ,d 3 k' 
or J 



dm k 
dr 



(14) 



Hence Eq. becomes 



dN k , 3N k 



dt y dr 



v g -^-+fj,k — / k'N k , d A k 



Or 



,\ dNk 
dk 



= 0, (15) 



which in the one-dimensional case reduces to 



dNk +v dNk 
dt 9 dx 



Hk J k'N k ' dk' 



,\ dNk 
dk 



0. (16) 



A similar equation may of course be derived for the in- 
tensity I k . 

Equation l|16(l gives the evolution of high-frequency 
photons on a slowly varying background radiation fluid, 
and it may be used to analyse the long term behaviour 
of amplitude modulated intense short incoherent laser 
pulses. 

We have considered the nonlinear propagation of ran- 
domly distributed intense short photon pulses in a pho- 
ton gas. The photon-photon interactions are described by 
means of a QED model, in which an ensemble of incoher- 
ent photon pulses is governed by a wave kinetic equation 
where the coupling between the intense photon pulses 
and the sound-like waves of the photon gas is due to slow 
variations of the sound wave energy distribution. The in- 
tense photon pressure, in turn, modifies the sound wave 
propagation. The wave kinetic and the sound wave equa- 
tions form a closed system, which has been used to derive 
a dispersion relation. By choosing appropriate spectra 
for short pulse photons, we analyze the dispersion rela- 
tion to show the existence of new classes of modulational 
instabilities. The latter can cause fragmentation of in- 
coherent photon pulses in astrophysical contexts and in 
forthcoming experiments using very intense short laser 
pulses. 
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